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Abstract 

Cosmic rays of energies larger than the Greisen-Zatsepin-Kuzmin 
(GZK) cutoff may be neutrinos if they acquire strong interactions due 
to a "precocious unification" of forces. A scenario for this to happen 
is outlined. There is no contradiction with precision measurements 
carried out at LEP and SLAC. Observable consequences at LHC and 
future neutrino detectors arc discussed. 



PACS 13.85.T, 13.15, 14.60.S 



A substantial number of cosmic ray events has been detected in which the 
primary energy appears to exceed the Greisen-Zatsepin-Kuzmin (GZK) cut- 
off, see e.g. ||, ||[ |[] and Szabelski's recent review ||. The onset of the 
GZK cutoff itself is somewhat uncertain: the energy of the rapid turnover of 
the spectrum depends on several details (for instance, on the injection spec- 
trum, etc. ), see ref. @ for a modern treatment. Nevertheless, it is unlikely 
that all the events reported can be explained by fluctuations in the shower 
development. 
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A number of explanations of the "anomalous" events has been offered; 
for an overview, cf. the proceedings of the 1997 University of Maryland 
workshop 0. In a recent article, however, Burdman, Halzen and Gandhi 
ref. || conclude that none of the explanations offered in the literature is a 
convincing one. 

Given this situation, we reexamine a proposal put forward some time 
ago, []|, |ip| . In those papers, we proposed that at sufficiently high energies, 
neutrinos (in general, leptons) acquire some unspecified strong interaction 
and cause (possibly) post-GZK showers. Since neutrinos are neutral, have 
small magnetic moments and at low energies they have nothing but weak 
interactions, they can freely propagate through the 2.7°K background. Con- 
sequently, they can reach us from cosmological distances. 

There were two weaknesses inherent in this suggestion. First, no mech- 
anism has been offered as to how the neutrinos get their strong interaction. 
Second, due to the fact that we used a sharp (0 function) threshold to turn 
on the strong interaction, the resulting amplitude violated unitarity, cf. ||. 

The scenario we present here is based on recent work designed to over- 



come the hierarchy problem. The authors of refs. pi] , 12, 13, 14, 15, |l£|] 
conjecture that unification may take place at a much lower energy than 
originally suspected, perhaps at a few TeV ( "precocious unification"). Even 
though there are some problems with this approach (e.g. the long lifetime 
of the proton is still lacking a convincing explanation), it offers some inter- 
esting possibilities regarding the nature of the physics beyond the Standard 
Model. 

The high unification masses obtained within the Standard Model and 



its minimal supersymmetric extensions, see ]17| , 18, 19, ^(J, is the result of 
the assumption of a "desert" between - approximately - the weak scale and 
the GUT scale. If one wants to accomplish unification at a lower energy, 
one has to postulate the existence of a rapidly increasing density of states 
somewhere above the weak scale. 

For the purpose of the present paper, it is irrelevant whether the rapid 
increase of the density of states is due to the existence of extra dimensions, 
due to "stringy effects" becoming relevant at lower energies or to something 
else. 

Let us now assume that we have a rapidly increasing level density, which 
has approximately the same form as in string theories. We choose: 

d{m n )^(^\ expf^j , (1) 



see for instance [21]. Here, m n stands for the mass of the resonance at 
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the n th level of a string model, tuq is the characteristic energy scale of the 
model. A numerical prefactor has been omitted: it is model dependent and 
it it does not play a significant role in what follows. Similarly, the exponent 
a is model dependent. We found, however that varying the exponent does 
not affect the results significantly; for the sake of definiteness, we settled 
with a = 2. 

We quoted an expression corresponding to the asymptotic form of the 
level density in a string model. Due to the fact that at this stage we want 
to understand the qualitative behavior of the neutrino cross section in the 
scenario just outlined, the present expression should be adequate. In most 
string models, the value of the exponent p is 1/2. Nevertheless, we wanted 
to explore more rapidly rising level densities as well; all the estimates were 
carried out with p = 1/2 and p = 1. 

In order to estimate the neutrino-quark cross section, we take into ac- 
count the contribution of the resonances to the imaginary part of the forward 
neutrino-quark amplitude. The invariant Breit-Wigner formula is used for 
a single resonance. In this way, we get the contribution of a single resonance 
of mass squared s n at level n of a string model: 

1 sT s l/2 
ImBn = fi^5 (2) 

In this equation s stands for the CM energy squared of the neutrino-quark 
system. 

The total cross section due to the "new physics" is then obtained by 
multiplying eq. (||) with the level density and summing over the levels. We 
assume that T n sl/ 2 grows linearly with s n . viz. T n sl/ 2 = jos n and that 
the resonances lie on a linear Regge trajectory, s n = sq(tl — 1). Further, 
we identify sq = tuq in eq. (|l|). (In a string model, both quantities are 
related to the string tension; there may be prefactors of order unity which we 
ignore here.) The linear growth of the width can be understood in intuitive 
terms by noticing that a resonance preferentially decays into channels with 
the largest phase space available. The number of such channels is roughly 

1/2 

proportional to the mass of the resonance, hence Y n Sn oc s n . (It is known, 
for instance that the total widths of baryon resonances grow approximately 
linearly with their masses cf. p2]].) 

We now notice that 70 is likely to be rather small, cf. (ll]]. Hence, the 
quantity ImB n may be replaced by a (5-function and the summation over 
the levels by an integration. In this approximation the neutrino-quark cross 
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section becomes: 

& fa — d (s) (3) 
so 

Finally, one has to integrate eq. (|3|) over the momentum distribution of 
quarks within the target nucleon, using the relation s = xs. Little or nothing 
is known about the evolution of structure functions into a region of rapidly 
increasing level density. In order to get some idea, we chose a "generic" 
structure function, of the form: 

S(x) = Ax' 1 (1 - xf (4) 

Inspired by the Duke-Owens parametrization of the structure functions, 
cf. [g3[ we chose A = 2, 7 = 0.6, 5 = 3.5. None of the results turned out 
to be very sensitive to the precise choice of these parameters. The final 
expression of the neutrino-nucleon cross section is given by: 

a = / dxS(x)a(xs) (5) 

Jx 

The infrared cutoff was chosen as xq = sq/s; in this way, the "new physics" 
begins to manifest itself for s ^ sq. In order to test the scheme developed, 
let us assume that we want the cross section to grow to a fa 1/ A.q CD around 
the GZK cutoff, say, s = 2 x 10 4 TeV 2 . (This corresponds to a laboratory 
energy, E L fa 10 19 eV.) On taking A QC d = 200MeV, one gets s ~ 3TeV for 
p = 1/2 in eq. (JlJ) and so ~ 30TeV for p = 1, respectively. In the following 
figures we display the cross sections resulting from eq. (p) taking p = 1/2 
and p = 1 in eq. ([!]). 

One sees that, for all practical purposes, the neutrino-nucleon cross sec- 
tion is dominated by the exponential growth of the level density; hence, the 
choice of the exponents in eq. (|l]) and eq. (||) is not critical. 

The cross sections rise very rapidly and, as it is always the case with 
tree amplitudes containing high spin particles in intermediate states, it will 
violate unitarity at a certain energy. The value of that energy is, however, 
difficult to determine at this stage. This is due to the fact that, like in any 
string model, at level n, resonances with spins < s < n are exchanged. 
(In some models the lower limit may be 1 and the upper limit nil, but 
this fact does not affect the general situation.) As a consequence, the usual 
semiclassical estimates of the unitarity limits of cross sections are not directly 
applicable. Only by calculating loop corrections will one be able to estimate 
the unitarity limits. 
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Figure 2: Same as Fig. 1, but with p = 1, s Q ' = 30TeV 
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It is to be emphasized that individual partial waves need not violate 
the unitarity bounds. In order to illustrate this point, consider a simple 
Veneziano amplitude, 

A(s,t) = B(l - a(s),l - a(t)) , (6) 

where B(x,y) is the Euler beta function and a is a linear Regge trajectory. 
This amplitude has a level density asymptotically described by eq. ([l]) with 
p = 1/2. By using Stirling's formula, one obtains an asymptotic estimate 
for the partial wave amplitudes. In particular, one finds for the S wave 
amplitude: 

Ao(s) ~ (2vra( S ))- 1/2 (a(s) » 1) . (7) 

Clearly, the bound given by |exp(i5o) sin6o| < 1 is satisfied asymptotically. 
(In order to arrive at the result given by eq. ([?]), one has to avoid the poles 
of the r function, e.g. by sending s to infinity along a ray in the complex 
plane, s = \s\ expicp, <p ^ 0.) 

Are the precision results of the Standard Model affected? This is an 
important question: the effects of the "new physics" are expected to be ob- 
servable even below the characteristic energy scale, cf. Goldberg and Weiler, 
ref. p3]. As in discussing the unitary bounds before, we are unable to es- 
timate loop effects at this stage. However, one can estimate final state 
interactions due to the "new physics" . We use the scattering length approx- 
imation to the I th partial wave, 

£^ +1 )cot^ - (8) 

This should be reasonably accurate: typical CMS energies at precision mea- 
surements carried out at LEP and SLAC are of the order of lOOGeV, whereas 
the characteristic energy of the "new physics" is on the TeV scale. A rea- 
sonable estimate for the scattering length is 1/a ~ \/sp- The effect of the 



final state interactions is given by the formula, cf. [25]: 



wi 



In eq. (^), w is a transition probability in a given partial wave I, (for instance, 
Z — ► vv) and w° is the same transition probability calculated ignoring the 
final state interaction. On using the estimate */sq ~ 3TeV, one finds that 
with k ~ 45GeV (half of the mass of the Z), the S wave final state "en- 
hancement factor" differs from unity by about 10 -4 or so. This is to be 
compared with a typical error of a precision measurement (for instance, the 
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decay width of a weak gauge boson) which is of the order of 0.1%, [26]. 
The effect is, of course, much smaller in higher partial waves or for a higher 
characteristic energy, cf. eq. (0). 

The scenario outlined here has some observable consequences. They, in 
turn, depend on whether the characteristic energy is of a few TeV or a few 
tens of TeV. For a characteristic energy of a a few TeV, 

• there should be measurable deviations from standard model predic- 
tions in the decay rates of weak gauge bosons if the accuracy of the 
measurements can be increased by about an order of magnitude, 

• one expects spectacular phenomena of new particle production and/or 
strong violation of Feynman scaling due to the rapid rise of the level 
density at LHC. Due to the precocious unification, one expects a co- 
pious production of leptons as well as hadrons. 

For a high characteristic energy (~ 30TeV), this phenomenon may not take 
place at the LHC. However, new phenomena will be observed at nonaccel- 
erator experiments, such as neutrino telescopes and at orbiting detectors 
(OWL, Airwatch). 

• It was pointed out that orbiting detectors should be sensitive to neu- 
trino interactions, see [27]. If the cross section of neutrino interactions 



shows a rapid rise, there should be a corresponding rise of the im- 
pact parameter of the incident neutrino with respect to the center of 
the Earth at which showers generated by them can be observed. Ac- 
cordingly, there should be a cutoff in the spectrum of upward going 
neutrinos observed in underground, under water or under ice detec- 
tors. The neutrino induced showers develop rapidly within the Earth 
and they degenerate by the time they would reach the detector. By 
contrast, the number of showers of grazing incidence should increase. 

This scenario resolves the "energy crisis" caused by some schemes pur- 
porting to explain the highest energy cosmic rays, such as ||. It 
is difficult enough to construct mechanisms by means of which pro- 
tons are accelerated to energies of the order of 10 20 eV, cf. Norman 
et. al. , ]^| . If one prefers the highest energy cosmic rays to be neu- 
trinos originating from the decay of pions and kaons, one needs protons 
to be accelerated to energies a few orders of magnitude even higher: at 
such energies, several hundred light hadrons are produced and, on the 
average, the available primary energy is shared equally between them. 
By contrast, in a scenario involving precocious unification, once the 
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CMS energy of the accelerated protons in an active galactic nucleus 
or in a similar site of intense proton acceleration reaches the charac- 
teristic energy, neutrinos are produced at multiplicities comparable to 
hadrons. Hence, there is no need to postulate proton energies several 
orders of magnitude larger. 



We thank L. Madansky and C. Norman for several conversations on the 
highest energy cosmic rays. 
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